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Abstract 

In this paper, we study the effect of feedback on two-user MIMO interference channels. The capacity region of 
MIMO interference channels with feedback is characterized within a constant number of bits, where this constant is 
independent of the channel matrices. Further, it is shown that the capacity region of a MIMO interference channel 
with feedback and its reciprocal interference channel are within a constant number of bits. Finally, the generalized 
degrees of freedom region for the MIMO interference channel with feedback is characterized. 
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I. Introduction 

Wireless networks with multiple users are interference-limited rather than noise-limited. Interference 
channel (IC) is a good starting point for understanding the performance limits of the interference limited 
communications [1-7]. Feedback can be employed in ICs to achieve an improvement in data rates [8-13]. 
However, most of the existing works on ICs with feedback are limited to discrete memoryless channels or to 
single-input single-output (SISO) channels. This paper analyzes multiple-input multiple-output (MIMO) 
Gaussian interference channels with feedback. In particular, we focus on a two-user MIMO IC with 
feedback. 

The authors of [8] considered a SISO Gaussian IC with feedback and found its capacity region within 
two bits. It was shown that the capacity regions of Gaussian ICs increase unboundedly with feedback 
unlike the Gaussian multiple-access channel where the gains are bounded [14]. The degrees of freedom 
for a symmetric SISO Gaussian IC with feedback is also found in [8]. In this paper, we find an outer 
bound and an inner bound for the capacity region that different by a constant number of bits, and also 
evaluate the generalized degrees of freedom (GDoF) region for a general MIMO IC with feedback. 

The first main result of the paper is the characterization of the capacity region of a MIMO IC with 
feedback within JVi + N 2 + max(7Vi, N 2 ) bits, where Ni and N 2 are the numbers of receive antennas 
at the two receivers. An outer-bound is obtained by first outer bounding the covariance matrices of 
both input signals and representing the outer bound as a region in terms of the cross-covariance matrix 
between the two input signals. This is further outer-bounded by a larger region that does not involve 
the cross-covariance matrix. The achievability strategy is based on block Markov encoding, backward 
decoding, and Han-Kobayashi message-splitting. This achievable rate and the outer bound are within 
Ni + N 2 + max(7Vi, N 2 ) bits of each other thus characterizing the capacity region of the two-user IC 
within constant number of bits where the constant is independent of the channel matrices. The achievability 
scheme that is used to prove the constant gap result assumes that the transmitted signals from the two 
transmitters in a time-slot are uncorrelated, unlike [8] where the signals were assumed correlated in the 
achievability. Thus, our achievable rate region is within 3 bits rather than 2 bits as in [8] of the capacity 
region of a SISO IC with feedback. An achievability scheme without correlated inputs was also shown 
to achieve within constant gap of the capacity region in [12] for a SISO IC with feedback. However, our 
gap between the inner and the outer bounds is smaller as compared to [12]. 

We note that the achievability strategies for a SISO IC in [8, 12] emphasize that the private part from 
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a transmitter using the Han-Kobayashi message splitting is such that it is received at the other receiver 
at the noise floor. However for a MIMO IC with feedback, it is not clear what its counterpart would be. 
The Han-Kobayashi message splitting used in this paper gives the notion of receiving the signal at the 
noise floor for a MIMO IC with feedback. Many matrix based results are derived in this paper to show a 
constant gap between the outer and the inner bounds of the capacity region of a MIMO IC with feedback, 
which may be of independent interest. 

The second main result of the paper is to show that the capacity region of a MIMO IC with feedback 
and that of its corresponding reciprocal channel are within constant number of bits of each other, where 
the constant is independent of channel matrices. The reciprocal IC was considered in [4], where the 
authors showed that the capacity region of a MIMO IC without feedback is within constant number of 
bits of its corresponding reciprocal IC. This paper shows that the constant gap between a MIMO IC and 
its reciprocal channel also holds in the presence of feedback. 

The third main result of the paper is a complete characterization of the GDoF region of a general MIMO 
IC with feedback when the average signal quality of each link, say p {j for link from transmitter % to receiver 
j, varies with a base signal-to-noise ratio (SNR) parameter, say SNR, as limsMR^oo = <*ij, where 

aij can be different for each link with i,j e {1,2}. In other words, the average link quality of each link 
can potentially have different exponents of a base SNR. As a special case, we consider a symmetric IC 
where the number of antennas at both transmitters is the same, the number of antennas at both receivers 
is the same, and the SNRs for the direct links and the cross links are SNR and SNR", a > 0, respectively. 
We find the GDoF (the maximum symmetric point in the GDoF region) for a given a and show that the 
GDoF is a "V"-curve rather than a "W"-curve corresponding to the GDoF without feedback as in [5]. 
Similar result was obtained for a SISO IC in [8] while this paper extends it to a MIMO system. 

The reminder of the paper is organized as follows. Section II introduces the model for a MIMO IC 
model with feedback, reciprocal IC and the GDoF region. Sections III and IV describe our results on the 
capacity region and the GDoF region respectively. Section V concludes the paper. The detailed proofs of 
various results are given in Appendices A-E. 

II. Channel Model and Preliminaries 

In this section, we describe the channel model considered in this paper. A two-user MIMO IC consists 
of two transmitters and two receivers. Transmitter % is labeled as Tj and receiver j is labeled as Dj for 
i,j G {1,2}. Further, we assume Tj has M { antennas and Dj has N { antennas, i e {1,2}. Henceforth, 
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such a MIMO IC will be referred to as the {M u N u M 2 , N 2 ) MIMO IC. We assume that the channel 
matrix between transmitter Tj and receiver Dj is denoted by G C NjXMi , for i,j G {1,2}. We shall 
consider a time-invariant or fixed channel where the channel matrices remain fixed for the entire duration 
of communication. We also incorporate a non-negative power attenuation factor, denoted as p i:) , for the 
signal transmitted from Tj to D^. At time-instant t, transmitter Tj chooses a vector X^t) G C MlXl and 
transmits y/P^X^t) over the channel, where Pj is the average transmit power at transmitter Tj, and 

tr(Qjj) < 1, (1) 

where Q { j = E(XjXj) for i,j G {1,2}. 

In this paper, we denote the conjugate transpose of the matrix S as Further, A<B means B — A is 
a positive semi-definite (p.s.d.) matrix and A >z B if B -< A. The identity matrix of size s x s is denoted 
by I s . Further, we define x + = max(x, 0). 

By definition of Q^, we see that Qij = Q^. We will sometimes denote Q = Qi 2 when it does not lead 
to confusion. 

In addition, since Qu is p.s.d., its eigenvalues are positive. tr(Qjj) is sum of eigenvalues and (1) implies 
that it is less than or equal to 1. So, all eigenvalues of Qu are less than or equal to 1. According to Theorem 
7.7.3. of [15] we get Qu ^ /. Moreover, we have ^ QijQlj d I, where ■< QijQlj results from the fact 
that every matrix in the form of AA^ in p.s.d. and QijQlj d I results from t?(QijQij^) = tr(Q ii )tr(Qjj) < 
1 which gives QijQ^ -< I with a similar argument as we had for Qu. 

The received signal at receiver Dj at time instant t is denoted as Y^t) for % e {1, 2}, and can be written 

as 

Yxit) = ^T 1 H 11 X 1 (t) + v ^ 2 - 1 H 21 X 2 (t) + Z 1 (t), (2) 
Y 2 (t) = y/p^H 12 X 1 (t) + y /p^H 23 X 2 (t) + Z 2 (t), (3) 

where Zi(t) e C NiXl is independent and identically distributed (i.i.d.) CN(0, ijvj (complex Gaussian 
noise), p u is the received SNR at receiver Dj and pij is the received interference-to-noise-ratio at receiver 
Dj for i,j G {1,2}, % 7^ j. A MIMO IC is fully described by three parameters. The first is the number 
of antennas at each transmitter and receiver, namely (Mi, N±, M 2 , N 2 ). The second is the set of channel 
gains, H = {H n , iJ 12 , H 21 , H 22 }. The third is the set of average link qualities of all the channels, 
P = {pn, P12, p 2 i, P22}- We assume that these parameters are known to all transmitters and receivers. 
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For MIMO IC with feedback, the transmitted signal at Tj, Xi(t), is a function of the message Wi 
and the previous channel outputs at receiver Dj for % e {1, 2}. Thus, the encoding functions of the two 
transmitters are given as 



where f it is the encoding function of Tj, Wi is the message of Tj and F/ -1 = (F(l), Yi(t — 1)). 
Let us assume that Tj transmits information at a rate of R { to receiver Dj using the codebook C i:Tl of 
length-n codewords with \C iyn \ = 2 nRi . Given a message m 8 e {1, . . . , 2 ni?i }, the corresponding codeword 
X"(mj) e C ijn satisfies the power constraint mentioned before. From the received signal Y™, the receiver 
obtains an estimate ml of the transmitted message m 8 using a decoding function. Let the average probability 
of error be denoted by e iy1l = Pr( ^ mi). 

A rate pair (R 1: R 2 ) is achievable if there exists a family of codebooks C^ n , i = {1, 2} n and decoding 
functions such that maxj{ej jn } goes to zero as the block length n goes to infinity. The capacity region 
C(H,p) of the IC with parameters H and p is defined as the closure of the set of all achievable rate 
pairs. 

Consider a two-dimensional rate region C. Then, the region C © ([0, a] x [0,6]) denotes the region 
formed by {(Ri, R2) '■ Ri, R2 > 0, ((i?i — a) + , (i? 2 — b) + ) e C} for some a, 6 > 0. Similarly, the region 
C ([0, a] x [0, b]) denotes the region formed by {(R u R 2 ) : R ± , R 2 > 0, ((i?i + a) + , (i? 2 + &) + ) G C} for 
some a, b > 0. Further, we define the notion of an achievable rate region that is within a constant number 
of bits of the capacity region as follows. 

Definition 1. An achievable rate region A(H, p) is said to be within b bits of the capacity region if 



In this paper, we will use the GDoF region to characterize the capacity region of the MIMO IC with 
feedback in the limit of high SNR. This notion generalizes the conventional degrees of freedom (DoF) 
region metric by additionally emphasizing the signal level as a signaling dimension. It characterizes the 
simultaneously accessible fractions of spatial and signal-level dimensions (per channel use) by the two 
users when all the average channel coefficients vary as exponents of a nominal SNR parameter. Thus, we 
assume that 



Xi(t) 



MW.X' 1 ), ie{l,2}, 



(4) 



A(H,p) QC(H,p) and A(H,p)®([0,b]®[0,b])DC(H,p). 



log(SNR)-^oo log(SNR) 



lim 



(5) 
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where e M + for all i,j G {1, 2}. In the limit of high SNR, the capacity region diverges. 

The GDoF region is defined as the region formed by the set of all (d 1 ,d 2 ) such that (c? 1 log(SNR) — 
o(log(SNR)), d 2 log(SNR) — o(log(SNR))) 1 is inside the capacity region. Thus, the GDoF is a function 
of link quality scaling exponents a^. We note that since the channel matrices are of full ranks with 
probability 1, we will have the GDoF with probability 1 over the randomness of channel matrices. 

The property of maintaining the same performance even if the direction of information flow is reversed 
is known as the reciprocity of the channel. For a MIMO IC with parameters (Mi, N 1: M 2 , N 2 ), H = 
{iJn, #12, #21, #22}, and p = {pn, P12, p 21 , p 22 }, the reciprocal MIMO IC has parameters (N 1 , Mi, N 2 , M 2 ), 
H R = {Hf 1: H 21 ,Hf 2 , # 2 T 2 }, and p R = {p n , p 21 , p l2) p 22 }. 

III. Capacity Region of MIMO Interference Channel with Feedback 

In this section, we will describe our results on the capacity region of the two-user MIMO IC with 
feedback. 

Our first result gives an outer bound on the capacity region of the two-user MIMO IC with feedback. 
Let 1Z (Q) be the region formed by (Ri,R 2 ) satisfying the following constraints for some covariance 
matrix Q = E[XiXj]: 

'a = o(log(SNR)) indicates that lim S NR->cx> wsnr) = °- 



R-y < \0gdet(I Nl +pnHuHl 1 +p2iH 21 Hl 1 + ^^H n QHl 1 + ^^ 1 H 21 Q j 'Hl 1 ), (6) 
R 2 < log det(/AT 2 + p 22 H 22 H\ 2 + p 12 H 12 H\ 2 + ^ p 22 p\ 2 E 22 Q x E\ 2 + ^/ p 22 pi 2 H 12 QH\ 2 ), (7) 
R-y < log det (l N2 + p 12 H 12 H\ 2 - p 12 H 12 QQ x + log det ^I Nl + puHuH^ - 

In 2 + Pi2H\ 2 H\ 2 ^/p~nHi 2 Q 



y/pnPi2HuH\ 2 y/p^HnQ 



12 



l M 2 



(8) 



/»'•> < 



< log det (l Nl + p 2 iH 2l H ] 21 - p 2 iH 2l Q ] QH ] 2 ^j + log det I Ni + p 22 H 22 H ] 22 - 



+ P2lH 21 Hl 1 sfp^ll-iyQ} 
y/P~2lQH 2 i I Mi 



-1 



V 'P22P2\H 2 i H\ 2 
y/p~22~QHl 2 

Ry + R 2 < log det (l N2 + p 22 H 22 H\ 2 + p l2 H 12 H{ 2 + ^ 'p 22 p 12 H 22 Q x H{ 2 + y/p 22 pi 2 H 12 QH< 
+ log det ^ + pnHyyH^ - ^/piJp^HnH^ v /pTT#nQ ] 



(9) 



^jv 2 + Pi2Hi 2 H\ 2 y/pi2~H 12 Q 
\fpv2Q ] H\ 2 I M2 



-, -l 



y/PnpuH^H^ 



(10) 



i?i + i? 2 < log det (/ M + pnHyyHl + p 21 H 2l H ] 21 + y/p^H^QH^ + y/p^H 21 Q ] Hl^j 
+ log det Jjv 2 + p 22 H 22 H\ 2 - ^fp^2~p2~xH 22 H\ x ^/p^H 22 Q^ 



Ini + P2\H 21 H 21 Jp~2\H 21 Q^ 



-. -l 



21 



'Mi 



V 'p22p2lH 2 \ H 22 

Jp~*iQH 22 



(ID 



Further, let 72. be the convex hull of 1Z (Q) for all covariance matrices Q. The following theorem outer 
bounds the capacity region of the two-user MIMO IC with feedback. 



Theorem 1. The capacity region of the two-user MIMO IC with perfect feedback Cfb is bounded from 
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above as follows 

C FB C TZ . (12) 



Proof: The proof is given in Appendix A. ■ 
From the definition of 1Z (Q), by substituting Q = and after some simplifications, we get that 1Z (0) 
is the region formed by i? 2 ) satisfying the following 

i*i < logdet(/ JVl + Pl iH ll Hl l + p 2l H 21 Hl), (13) 
i? 2 < log det(7iv 2 + p 2 2^22 Hl 2 + p 12 H 12 Hl 2 ), (14) 
i?i < log det (7jv 2 + Pi2#i2#i 2 ) + log det(7 JVl + p\\E\\E[ x - 

^l^pT 2 H ll H\ 2 (I N2 + p 12 H l2 H\ 2 y l ■y/p^H^H^), (15) 
i? 2 < log det + p 2 iH 2 iH\^j + log det(ijv 2 + p 22 H 22 H\ 2 - 

\fp22p2iH 22 H\ l (I Nl + p 2 \H 2X H\^ x p 22 p 2X E 2l H\ 2 ), (16) 
i?i + i? 2 < log det (/jv 2 + p 22 H 22 H\ 2 + p 12 H 12 H\ 2 ^j + log det(/ JVl + PuH^H^ - 

y/pl^HnHl^I^ + p l2 H l2 H\ 2 )-' VpHpT^tfn), (17) 
i?! + i? 2 < log det (7^ + puHnHl + p2iH 2l H\^j + log det(/jv 2 + p 22 H 22 H\ 2 - 

^/p22p2\H 22 Hl l (I Nl + p 2 \H 2X H\^ x p 22 p 2X E 2l H\ 2 ). (18) 

The following result gives an inner bound to the capacity region of the two-user MIMO IC with 
feedback. 

Theorem 2. 77ze capacity region for the two-user MIMO IC with perfect feedback Cfb is bounded from 
below as 

C FB 2 1l o (0) Q ([0,7V! + 7V 2 ] x [0,^ + N 2 \). (19) 



Proof: The proof is provided in Appendix B. ■ 
The inner bound uses the achievable region for a two-user discrete memoryless IC with feedback as in 
[8]. The achievability scheme employs block Markov encoding, backward decoding, and Han-Kobayashi 
message- splitting. This result for a discrete memoryless channel is extended to MIMO IC with feedback 
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using a specific message splitting by power allocation. The transmitted signal from Tj, X it is given as 

Xi = Xi P + Xi U , (20) 

where X ip indicate the private message of Tj, and X iu is the public message of Tj. We assume that (X ip 
and X iu pairs are independent for i = 1,2.. However, these transmitted signals are correlated over time 
due to block Markov encoding. The private signal X ip is chosen to be X ip ~ CN (0, K Xip ), and the public 
signal X iu is chosen to be X iu ~ CN (0, K x . u ), where 

^x ip = /m, - x'Pull^lx, + PijHijH^^pTjHij, (21) 

and 

K Xiu =I Mi -K Xip , (22) 

for i G {1,2}. 

We will show in Appendix B that the power allocation is feasible by showing K Xip >z and K Xiu >z 0. 
Further, this message split is such that the private signal is received at the other receiver with power 
bounded by a constant. More specifically, we have p,j i Hj i K Xtp H^ i -< I N . thus showing that the effective 
received signal covariance matrix at receiver Dj corresponding to the private signal from transmitter Tj 
is at the noise floor. 

This power allocation is different from that given in [8] even for a SISO channel. Note that the power 
split levels in the achievability scheme of [8] do not sum to 1 and thus do not satisfy the total power 
constraint. For the special case of SISO IC with feedback, the above gives a fix to the results in [8]. This 
power allocation assumes uncorrelated signals transmitted by the two users at each time-slot. The authors 
of [12] also used uncorrelated signals for SISO but had a larger gap between the inner and outer bounds 
for SISO IC with feedback than that achieved by our achievability strategy. 

Having considered the inner and outer bounds for the capacity region of the two-user IC with feedback, 
the next result shows that the inner bound and the outer bound are within Ni + N 2 + max(A r 1 , N 2 ) bits 
thus finding the capacity region of the two-user IC with feedback, approximately. 

Theorem 3. The capacity region for the two-user MIMO IC with perfect feedback Cfb is bounded from 
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above and below as 



K o (0) © ([0, Ni + N 2 ] x [0, N x + N 2 }) C C FB C ft o (0) © ([0, JVJ x [0, N 2 ]) 



(23) 



where the inner and outer bounds are within N\ + N 2 + max (JVi, iV 2 ) tes. 

Proof: The inner bound follows from Theorem 2. For outer bound, we outer-bound the region 1Z (Q) 
as n o {Q) C 7e o (0)©([0, iVi] x [0, JV 2 ]) in Appendix C. Hence, ft C ft o (0)©([0, iVi] x [0, N 2 \). Thus, using 
Q = in 1Z (Q) gives an approximate capacity region with the approximation gap as in the statement of 



The authors of [8] found the capacity region for the SISO IC with feedback within 2 bits. The above 
Theorems generalizes the result to find the capacity region of MIMO IC with feedback within iVi + iV2 + 
max(iVi, N 2 ) bits. Note that the approximate capacity region without feedback in [4] involves bounds on 
2Ri + R 2 which do not appear in our approximate capacity region with feedback. In addition, in [8], the 
approximate capacity region for the SISO IC with feedback involves the cross-covariance matrix of the 
inputs in the inner and outer bounds, whereas our approximate capacity region for the MIMO IC with 
feedback does not involve cross-covariance matrix of the inputs. 

Figure 1 gives a pictorial representation for the result of Theorem 3. The inner and the outer bounds 
for the capacity region for MIMO IC with feedback are within a constant number of bits from the region 
TZ o (0) and thus the inner and outer bound regions are within a constant number of bits of each other. 



Fig. 1. Inner and outer bounds for the capacity region of MIMO IC with feedback are within a constant number of bits. The arrows from 
the corners A and B in TZ o (0) toward their respective corners on outer bound have vertical length of Nj and horizontal length of 7V2. 
The arrows from the corners A and B in 1Z (0) toward their respective corners on inner bound have the vertical and horizontal length of 
Ni + N 2 each. 



the Theorem. 



R 




N +N N 

1 2 2 



Having characterized the approximate capacity region for the MIMO IC with feedback, we next explore 
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the relation of capacity region of the MIMO IC with feedback with that of the corresponding reciprocal 
MIMO IC with feedback. The next theorem shows that the capacity region of the MIMO IC with feedback 
is approximately the same as that of its corresponding reciprocal channel with feedback. 

Theorem 4. The capacity region for the two-user MIMO IC with feedback Cfb and the capacity region 
for its corresponding reciprocal IC Cp B with feedback are within constant gaps from each other. More 
precisely, the following expressions holds: 

n o (0)e([0,N 1 + N 2 ] x [0,N 1 + N 2 ]) CC FB C ft o (0)©([0,7Vi] x [0,iV 2 ]), (24) 
U o {0) © ([0, M 1 + M 2 ] x [0, M 1 + M 2 \) C C% C K o (0) © ([0, x [0, M 2 ]). (25) 



Then, we get 



C% © ([0, JVi + N 2 + Mi] x [0, Ni + N 2 + M 2 \) C C FB C 
C^ B © ([0, Mi + M 2 + iVi] x [0, Mi + M 2 + jV 2 ]), (26) 
C FB © ([0, Mi + M 2 + jVi] x [0, M 1 + M 2 + N 2 \) C C 
C FB © ([0, iVi + iV 2 + Mi] x [0, Ni + N 2 + M 2 ]). (27) 



Proof: In Appendix D, we show that the region TZ o (0) for the MIMO IC is the same as the 
corresponding region TZ^(0) for the corresponding reciprocal MIMO IC. Thus, (24)-(25) follow from 
Theorem 3. Moreover, (26)-(27) follow by simple manipulations on (24)-(25). ■ 
Thus, we see that the capacity region of a two-user MIMO IC with feedback and the corresponding 
reciprocal channel with feedback are within N x + N 2 + M x + M 2 + max (iV\ + M 1 ,N 2 + M 2 ) bits. 

IV. Generalized Degrees of Freedom Region of MIMO Interference Channel with 

Feedback 

This section describes our results on the GDoF region of the two-user MIMO IC with feedback. The 
GDoF gives the high SNR characterization of the capacity region. Since the inner and outer-bounds on 
the capacity region are within a constant gap, we characterize the exact GDoF region of the MIMO IC 
with feedback. 
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Define 

{min (u, Ui)a{ + min {{u — Ui) + , u 2 )a 2 , if ai > a 2 
(28) 
min (u, 112)0,2 + min ((u — u 2 ) + , u±)af : otherwise 

The following result characterizes the GDoF for general MIMO IC with feedback for general power 
scaling parameters o;^-. 

Theorem 5. The GDoF region of the two-user MIMO IC with feedback is given by the set of (d\,d 2 ) 
satisfying: 

audi < f(N 1 ,(a 11 ,M 1 ),(a 21 ,M 2 )), (29) 
a 22 d 2 < /(JV 2 ,(a22,M 2 ),(ai2,Mi)), (30) 
audi < ai 2 min (Mi,N 2 ) + anmin ((Mi - N 2 ) + , JVi) + 

(«n-«i 2 ) + (min(M 1 ,iV 1 ) - min ((M x - N 2 ) + , N x ) ), (31) 
a 22 d 2 < a 2 imin (M 2 , A^i) + a 22 min ((M 2 - iVi) + , iV 2 ) + 

(a 22 - a 21 ) + (min (M 2 , N 2 ) - min ((M 2 - N ± ) + , N 2 ) ) , (32) 
«iidi + a 22 d 2 < /(A^ 2 ,(a 22 ,M 2 ),(a 12 ,M 1 )) + Q;iimin((M 1 -^ 2 ) + ,^ 1 ) + 

(an - a 12 ) + (min (M 1? - min ((M, - N 2 ) + , TV,) ) , (33) 
ct xl d x + ct 22 d 2 < f(N 1 ,(a 11 ,M 1 ),(a 21 ,M 2 )) + a 22 mm((M 2 -N 1 ) + ,N 2 ) + 

(a 22 - a 21 ) + (min (M 2 , N 2 ) - min ((M 2 - N t ) + , N 2 ) ) . (34) 

Proof: According to Theorem 3, we can see that GDoF is equal to limsNR^oo 7£ o (0)/ logSNR. 
Appendix E evaluates the limit of 7?. o (0)/logSNR as SNR — > 00 to get the result as in the statement of 
the Theorem. ■ 
Since the capacity region of the MIMO IC with feedback and the corresponding reciprocal IC with 
feedback are within constant gap, the GDoF region of the MIMO IC with feedback and that of the 
corresponding reciprocal IC with feedback are the same, as given in the next corollary. 

Corollary 1. The GDoF region for the reciprocal IC with feedback is given by the set of (d±, d 2 ) satisfying 
(29>(34). 

We will now consider a special case of Theorem 5 where M x = M 2 = M, N x = N 2 = N, a n = a 22 = 
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1, and CK12 = CC21 = a. This MIMO IC is called a symmetric MIMO IC. We also define GDoF d as the 
supremum over all di such that (di, di) is in the GDoF region. The GDoF for the symmetric MIMO IC 
with feedback is given as follows. 

Corollary 2. The GDoF for a two-user symmetric MIMO IC with feedback for N < M is given as 
follows: 

GDoFp F — / N ~^N-Mr, <f«<l, (35) 
{ N(^)-t(2N-M)+, ifa> 1. 

SzVzce ?/ze expressions are symmetric in N and M by Corollary 1, the GDoF for M < N follow by 
interchanging the roles of M and N. 

Proof: For the symmetric MIMO IC, we have 

/(AU^M^feM,)) = f(N, (1, M) , (a, M)) 

= max(l, a) min(M, N) + min(l, a) min((iV - M) + ,M). (36) 

Let the GDoF be d. We will split the proof for N < M in two cases. 
Case 1 - a < 1: We will go over all equations (29)-(34) and evaluate them for the symmetric case with 
a < 1. Equations (29) and (30) can be simplified using (36) as follows 

d < max(l,a)min(M, iV) + min(l, a) min((iV - M) + ,M) 

= N. (37) 

Equations (31) and (32) can be simplified as 

d < amin(M,A0 + mm ((M - N) + , N) + (1 - a) + (min (M, N) - mm ((M - N) + , N) ) 

= aN + min((M - N), N) + (1 - a)iV - (1 - a)min ((M -N),N) 

= N + amin((M - N),N) 

= N + a(N- (2N-M) + ). (38) 
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Equations (33) and (34) can be simplified as 

d < -(max(l,a)min(M,iV) + min(l, a) min((iV - M) + , M) + min ((M - N) + , N) + 
2 

(1 -a) + (mm(M,iV) — min ((M — iV) + , N) )) 
= ^(7V + (1 — a)N + a min((M — N), N)) 
= N - -a(N - (N - (2N - M)+)) 

= N — —((2N — M) + ). (39) 

We note that the minimum of the right hand sides of (37), (38), and (39) would give us the GDoF. The 
minimum of these three terms is (39) which proves the result for a < 1. 
Case 2 - a > 1: 

In this case, equations (29) and (30) can be simplified as 

d < max(l,a)min(Af, JV) + min(l, a) min((iV - M) + ,M) 

= aN. (40) 

Equations (31) and (32) can be simplified as 

d < amin (M, N) + min ((M-N) + ,N) + 

(1 - a) + (min (M, N) - min ((M - A^) + , N) ) 
= aN + min^M- N),N). (41) 

Equations (33) and (34) can be simplified as 

d < i(max(l,a)min(M,A^) + min(l,a)min((A^-M) + ,M) + min((M-A^) + ,A^) + 
(1 -a) + (min(M,iV) — min ((M — iV) + , N) )) 
= ^(aN + (N - (2N - M) + )) 

= jV (« + l) _ h 2 N-M) + . (42) 
2 2 

We note that the minimum of the right hand sides of (40), (41), and (42) would give us the GDoF. The 
minimum of these three terms is (42) which proves the result for a > 1. ■ 
The authors of [5] found the GDoF for the two-user symmetric MIMO IC without feedback as follows 
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for iV < M (We can interchange the roles of N and M if iV > M.) 



GDoF NF 



N-a{2N-M) + , 
N — (1 — a)(2N — M) + , 
N- f(2iV-M) + , 
min{iV, iV(^) - i(2JV - M) + }, 



if < a < §, 
if \ < a < I 
if | < ol < 1, 
if a > 1 



(43) 



We note that the GDoF with and without feedback are the same for | < a < 1. Figure 2 compares 
the GDoF for the two-user symmetric MIMO IC with and without feedback. In Figure 2(a), the "W"- 
curve obtained without feedback delineates the very weak (0 < a < |), weak (\ < a < |), moderate 
(| < a < 1) 5 strong (1 < a < 3 — ^) and very strong (3 — ^ < a) interference regimes. In the presence 
of feedback, the "W"-curve improves to a "V"-curve which delineates the weak (0 < a < 1) and strong 
(1 < a) interference regimes for all choices of N and M. For ^ < N < M, we see that the GDoF 
with feedback is strictly greater than that without feedback for < a < 2/3 and for a > 3 — M/N. For 
N < M/2, we see that the GDoF with feedback is strictly greater than that without feedback for a > 2. 
The GDoF improvement indicates an unbounded gap in the corresponding capacity regions as the SNR 
goes to infinity. 

GDoF i 1 GDoF 



3N/2 




(a) f < N < M [a = 3 



(b) N < f 



Fig. 2. GDoF for symmetric MIMO IC with perfect feedback (PF), and no-feedback (NF) for (a) =f < TV < M, and (b) N < =f . 



Interestingly, from Figure 2(b) we can see that if we increase M when iV < 4^, the GDoF does not 
change. This can be interpreted as that while iV < =y, N act as a bottle-neck and increasing M does not 
increase the GDoF. As a special case consider a MISO IC for which we note that the GDoF is the same 
for all M > 2. Thus, increasing the transmit antennas beyond 2 does not increase the GDoF. However, 
increasing the transmit antennas from 1 to 2 gives a strict improvement in GDoF for all a > 0. Similar 
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result also holds for SIMO systems where increasing the receive antennas from 1 to 2 help increase GDoF 
while increasing the receive antennas beyond 2 does not increase the GDoF. 

V. Conclusions 

This paper gives the capacity region of a MIMO IC with feedback within Ni + N 2 + max(iVi, N 2 ) bits. 
The achievability is based on block Markov encoding, backward decoding, and Han-Kobayashi message- 
splitting. The capacity region for a MIMO IC with feedback is shown to be within a constant number of 
bits from the capacity region of the corresponding reciprocal IC. Further, the GDoF region for general 
MIMO IC is characterized. It is found that for symmetric IC, the GDoF form a "V"-curve rather than the 
"W"-curve without feedback. 

Appendix A 
Proof of Outer Bound for Theorem 1 

In this Appendix, we will show that Cfb Q Tlo{Q) for some covariance matrix Q = E[XiX|]. 

The set of upper bounds to the capacity region will be derived in two steps. First, the capacity region is 
outer-bounded by a region defined in terms of the differential entropy of the random variables associated 
with the signals. These outer-bounds use genie-aided information at the receivers. Second, we outer-bound 
this region to prove the outer-bound as described in the statement of Theorem 1. 

The following result outer-bounds the capacity region of two-user MIMO IC with feedback. 

Lemma 1. Let Si be defined as Si = ^pYjHijXi + Zj. Then, the capacity region of a two-user MIMO 
IC with feedback is outerbounded by the region formed by (Ri,R 2 ) satisfying 



Ri 


< 


KY 1 ) 


-Hz,), 


(44) 


R2 


< 


KY 2 ) 


-h(Z 2 ), 


(45) 


Ri 


< 


h(Y 2 


X 2 )-h(Z 2 ) + h(Y 1 \X 2 ,S 1 )-h(Z 1 ), 


(46) 


R2 


< 


h{Y 1 


X 1 )-h(Z 1 ) + h(Y 2 \X 2 ,S 1 )-h(Z 2 ), 


(47) 


R\ + R 2 


< 


h(Y 1 


S 1 ,X 2 )-h(Z 2 ) + h(Y 2 )-h(Z 1 ), 


(48) 


R\ + R 2 


< 


h(Y 2 


S 2 ,X 1 )-h(Z 1 ) + h(Y 1 )-h(Z 2 ). 


(49) 
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Proof: The proof follows on the same lines as the proof of Theorem 3 in [8], replacing SISO channel 
gains by MIMO channel gains and is thus omitted in this paper. ■ 

The rest of the section outer-bounds this region to get the outer bound in Theorem 1 . For this, we will 
introduce some useful Lemmas. 

The next result outer-bounds the entropies and the conditional entropies of two random variables by 
their corresponding Gaussian random variables. 

Lemma 2 ( [16]). Let X and Y be two random vectors, and let X G and Y G be Gaussian vectors with 
covariance matrices satisfying 



(50) 





X 




' X G ~ 


Cov 




= Cov 






Y 




Y G 



Then, we have 



h(Y) < h(Y G ), 
h(Y\X) < h (Y G J X G ) . 



(51) 
(52) 



The next result gives the determinant of a block matrix, which will be used extensively in the sequel. 

A B 



Lemma 3 ( [17]). For block matrix M = 



C D 



with matrices A, B, C, and D, we have: 



detM 



det A det(D - CA~ l B), if A is invertible, 
det D det (^4 - BD~ l C), if D is invertible. 



(53) 



Now, we introduce a Lemma that is a key result which will be used to upper-bound a conditional 
entropy term in this section and also to show an upper bound in Appendix C. 



Lemma 4. Let L(K, S) be defined as 



L{K,S) = K- KS{I N2 + S^KS) Vif, 



(54) 



for some Mi x M 1 p.s.d. Hermitian matrix K and some Mi x N 2 matrix S. Then if ■< K ± ■< K 2 for 
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some Hermitian matrices K\ and K 2 , we have 

L(K U S) ±L(K 2 ,S). (55) 

Proof: We note that since K is p.s.d., K + eI Ml is invertible for all e > 0. Given ^ K x -< K 2 , let 
F(e) = L(K 2 + eI Ml ,S) - L{K X + eI Ml ,S). We need to show that F(0) y 0. 

We first show that F(e) h for all e > 0. From Woodbury matrix identity (Appendix C.4.3 of [18]), 
we have that if A is invertible, (A + BDy 1 = A' 1 - A- l B(I + D A' 1 B)~ l D A' 1 . Thus, we have 
L(K + eI Ml ,S) = {{K + eI Ml )~ l + SS r y l by substituting A as (K + eI Ml )~\ B as S and D as 5 f in 
the above identity. 

Thus, F(e) = ((K 2 + e/ Ml ) _1 + SS^y 1 - ((K 1 + eI Ml )~ 1 + Stf*)" 1 . Since K x and X 2 are Hermitian 
p.s.d. matrices with K\ -< K 2 , it easily follows that F(e) y 0. 

Having shown that F(e) y for all e > 0, we will now prove the continuity of F(e) at e = 0. For this, 
we take the partial derivative of F(e) at e = and show that it is not unbounded thus proving that F(e) 
is continuous at e = 0. Thus, we have 

^ = ^L(K 2 + eI Ml ,S)-L(K 1 + eI Ml ,S) 
ae ae 

= ^L(K 2 + eI Ml ,S)-^L(K 1 + eI Ml ,S) (56) 
Thus, it is enough to show that lim e _> -^L(Ki + eI Ml ,S) is bounded. Thus, we have 

lim ^L{K t + eI Ml ,S) = lim ^-(K, + eI Ml - (K t + eI Ml )S(I Na + S\K t + e/ Ml )S)~ + eI Ml )) 

= I Ml - lim ^((Ki + eI Ml )S(I N2 + S\K t + eI M JS)~ 1 &(K i + eI Ml )) 
e-s>o ae 

= I Ml - S(I N2 + tfKi^&Ki - K t S(I N2 + S^KiSy 1 ^ 

+KiS(I N2 + S^KiSy^S^ + rfKiSy^Ki, (57) 

which is bounded. Thus, F(e) is continuous at e = 0. Further, since K x and K 2 are Hermitian, we see 
that F(e) is Hermitian and thus normal. From the Wielandt-Hoffman theorem [19], we note that the L 2 
norm of the difference in eigen-values (ordered in a particular way) of two normal matrices is bounded 
by the Frobenium norm of the difference of the two matrices. This shows that since F(e) y and 
F(e) — F(0) — > as e — > 0, we have that the eigen-values of F(e) approach the eigen-values of F(0) as 
e — > 0. Thus, all the eigen-values of F(0) are non-negative which proves that F(0) is positive semi-definite 
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thus proving the result. ■ 
The next three Lemmas outer-bounds entropy and conditional entropies of some random variables. 

Lemma 5. The entropy of the received signal at the i th receiver, h(Yi), is outer-bounded as follows 

h{Yi) < log det (l Ni + puHuHl + pyllyll], + x 'i>,,i> l ,ll„Q, ! H J , + y/PuPaHjiQ^H^ 

+ logdet(27re/ JVi ), (58) 

fori,] e {1,2}, i^j. 
Proof: 

KY,) < h(Yf) 

= logdet27reE[|y. G | 2 ] 

= log det 2vre (l N . + puH^Q^H^ + pyJl y,Q y s ll ], + x 'p,,pyll ,,Q,,U ], 

< log det 2ne (j Ni + puHuH^ + pyllyll + py.pyUy.Q.Jl ], + y/puPjiHjiQljHl^ 
= log det (j N . + PuHhHI + pyllyM], + x > papj.U nQ,;^ ,, + V 'PaPji H jiQlj H l) 

+ logdet(27re/jO, (59) 

where (a) follows from Lemma 2, and (b) follows from the fact that log det(.) is a monotonically increasing 
function on the cone of positive definite matrices and we have Qu ^ Im^m % for % e {1, 2}. 

Taking 27re out of the above determinant in the last part, gives the result as in the statement of the 
Lemma. ■ 

Lemma 6. The conditional entropy of the received signal at the i th receiver given the transmitted signal 
from the i th transmitter, hiY^Xj) is outer-bounded as follows 

h (Yi\Xi) < log det (j N . + PyMyll], - pjiHjiQljQijH}^ + log det 2ire (I Ni ) , (60) 

where Qij is the cross-covariance between Xi and Xj and Qa is the covariance matrix for Xj. 
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Proof: Let 



K a = E 



(61) 



where 



E P^ f ] = ^ + PuHuQuHI + pjiHjiQjjHjt + ^I'nl'jillnQijIIji + V 'PaPji H jiQ\j H iv (62) 



and 



K l2 4 E[X,X/] = g,,. 



(63) 



According to Lemma 2, we get 



MiTO < Mif |xf) 

= h{X?,Y?)-h{X?) 

= log det 2ire(Ka) — log det 27re(K i2 ) 

= logdet(^i) - logdet(K i2 ) + log det 2Ke(I Ni ). 



(64) 



Due to the reason that Q's elements are chosen from a continuous space, it is invertible with probability 
of one. In addition, according to Corollary 7.7.4(a) of [15], if we have Qu ^ Im^m^ Qu 1 h 1 m v 
Using Lemma 3 with M = K iX and A = K i2 , we get 



log det Kn 



(a) 



(b) 



log(det(E(X i X/)) det(E(r i F/) - E(r i X/)(E(X i X/))- 1 E(X i r/))) 
logdet(E(X,X i t )) + logdet(E(r i F. t ) - E(r,X i t )(E(X J X i t ))- 1 E(X,r/)) 
logdet(Qii) + logdet(/ iVl + p y M y.Q. Jl ], - pjiHjiQ^Q^QijHU 



< log det(Qii) + log det (7^ + p } ,\l y,\\ ], - ,) ji n i ;Q] j Q u Il i ji ). 



(65) 



where (a) is obtained by using (61) and some simplifications, and (6) follows from the fact that log det (.) 
is a monotonically increasing function on the cone of positive definite matrices and we have Qu ■< lu^Mi 
and Ql 1 h ^ Mi x Mi according to Corollary 7.7.4(a) of [15] for i, j e {1, 2} , % ^ j. 
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Substituting (65) in (64) gives the result as in the statement of the Lemma. ■ 

Lemma 7. The conditional entropy of the received signal at the i th receiver given Xj and S it h(Yi\Xj, Si) 
is outer-bounded as follows 

h (Yi | Xj, Si) < logdet ^I N . + PnHaHl - [ y /p^ j H ii H} j x 'p,, I !,,();, ] 



Inj + PijHijHjj ^/pijHijQi 

\/PijQij H ij 



-, -1 



y/PiiPij H ij H ii 



+ logdet 2ire (ijyj • 



(66) 



Proof: Let K i3 and K i4: be defined as follows 



and 



K 



i3 



E 



/PiiHaXi + Zi 

/~p~ijHijXi + Zj 
X, 



/PiiHaXi + Zi 

ffkjHijXi + Zj 
X, 



In, 

y/PiiPijHijQiiHii iNj + PijHijQiiHj- 
\/P~iiQijHii 



E 



fp~ijHijXi + Zj 
X, 



I PiiH'iiQij 
1 ' Pij HijQij 

Qjj 



PwHijXi + Zj 
X, 



iNj + PijHijQiiH\j \J Pij HijQij 



WjQljH, 



\jH\ 3 Q. 



33 



(67) 



(68) 
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Further, let Y{ = JpuHaXi + Zj. Then, 



h(Yi\Xj, Si) — h [^^fpiiHaXi + y/pjiHjiXj + Zi | Xj, ^JfJijHijXi + -Zj) 
— ^ {\Z~PuHiiXi + Zi | Xj, yfp^jHijXi + Zj) 



< WW,*?) 

= log det 27re(X i3 ) — log det 2ire(K i4 ) 

= log det (K a ) - log det (X i4 ) + log det 2ire(I Ni ), 



(69) 



where (a) follows from Lemma 2 by taking the two vectors Si and X, of lengths Nj and M, , respectively, 
together as a single vector of length of Nj + Mj and then, used Lemma 2. 
Substituting M = X i3 and D = K i4 in Lemma 3, we get 



log det (K i3 ) = log det (K i4 ) + log det y N . + puHuQuH^ 

~ ^l>iil>.jlliiQ::U,j y/piiHuQij ] [{K iA ) ] 1 

= log det (K u ) + log det I N . + puHuQaHt - 



\J PiiPijHij Qu Ha 
\/P~iiQijHii 



yf PiiPijHuQuH^j PiiHuQij 



-, -1 



\J PiiPijHij Qa 

VpTiQlM 



(70) 



+ PijHijQaHjj y/p~ijHijQij 

Note that since Qjj ^ /, using Lemma 3 we can see that Qjj = I outer-bounds the determinant of 

Inj + PijHijQuIljj y/p~ijHijQij 

y/PijQijHij Qjj 

Since B -< I implies ABA^ -< AA\ we have that Qjj = I outer-bounds the expression of the right 
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hand side of (70). Thus, 

log det (K i3 ) < log det (K i4 ) + log det ( I Ni + puHaQaHl - 



-, -1 



I X; + I'.jH ijQiill ;, y/f>ijHijQ 



^PiiQijHii 



WjQljHjj i Mj 

Next, we will show that Qu = I maximizes (71). 

Let us define S 4 ^if* , W 4 Q u - Q^-Q^ E A ^ + fi-t^-i and 

-l 



/ (S, Qu) — Qi 



QiiS Q 



QnS I Mi 



We can check that 



QljS i Mi 



E -ES^Q tj 
-QijSE I + Q^SE&Qa 



I 



Mi+Ni- 



Hence: 



/ (S, Qu) — Qi 



Qi 



QuS Q 



QiiS Qij 



E 



I Nj +&Q u S S^Qij 

QijS 

-ES^Q l3 

' ' ij 



S^Q l 

Ql 



^ PiiHaQij 

(71) 



-QljSE I + QljSE&Q 
= Qu ~ QiiSES^Qa + QaSES^QijQlj + QijQljSES^Qu — Q^Qlj — 

QijQij SES^QijQlj 

Qii QijQij (Qii QijQij)SES^ (yQa QijQij) 
— Qii — QijQij — (Qii — QijQ\j)S(I + S\Qu — QijQ\j)S) 1 S\Qu — QijQij) 

= w - ws(i + tfwsyi&w. 



(72) 



(73) 



(74) 



We know that W = Qa — QijQij d I — QijQij- So, according to Lemma 4 with K 1 as Qa — QijQij 
and K 2 as / — QijQij, we have / (S,Qu) ■< f(S,I Mi )- Thus, we use this outer-bound by replacing Q u 
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by I to get 



log det (K i3 ) - log det (K i4 ) 
< log det I N . + puHuHl - \ x /i>„ i > ij H,,H ;j x <p,,ll„Q IJ 



l\, + PijHijHij JWjHijQi 



-, -i 



WjQljHlj Im 3 



y/PiiPij H ij H ii 



Substituting this in (69), we get 



hiYilX^Si) < \ogdet(K l3 ) - log det (if*) + log det 2vre (I N .) 

< log det ^I Ni + Pii HiiHl- ^ l >„i>; l U„U ij ^i>;,U,,Q;, _ 

y/PiiPijHijHl 



Av: + pijH.jHjj yfpTjHijQij 
+ log det 2ne (I Ni ) . 



(75) 



(76) 



The rest of the section considers the 6 terms in Lemma 1 and outer-bounds each of them to get the 
terms in the outer-bound of Theorem 1 . 
First term: For the first term in Lemma 1, 

#i < M^l) - H z i) 

< log det (j Nl + pnHuH^ + p 21 H 21 H j 2l + a/PiiP 2 i#ii<2i2#21 + \ / Pnp2iH 2 iQ\ 2 Hl 1 ^ 

+ log det 2vre (I Nl ) - h{Z x ) 
( = } log det (l Nl + pnHnHl + p 2X H 2X E\ x + y/pl^HnQ^H^ + v / ) ^IT^# 21 Q t 12 # 1 t 1 ) , (77) 

where (a) follows from Lemma 5 and (b) follows from the fact that h(Zi) = logdet27re (ijvj. 
Second term: The second bound is similar to the first bound by replacing 1 and 2 in the indices. 
Third term: For the third bound in Lemma 1, it is sufficient to replace upper bounds of h (Y 2 \ X 2 ) and 
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h(Yi\X 2 , Si) from Lemma 6 and Lemma 7 as follows 



Ri < h(Y 2 \X 2 )-h(Z 2 ) + h(Y 1 \X 2 ,S 1 )-h(Z 1 ) 

(a) 



< log det (l N2 + p 12 H 12 H\ 2 - pi2H 12 Ql l Q2iH\ 2 ^ + log det 2ne (I N2 ) 



+ logdet I Nl +p 11 H u Hl 1 - 



V ' P\\PYiH\\tl\ 2 \fp~i\HiiQi2 



VpTiQ\ 2 hI 



In 2 + Pi2H\ 2 H l2 y/pi2~H 12 Q 12 

\fP~v2Q\2H\2 Im 2 
+ logdet27re (I Nl ) - h(Z 1 ) - h(Z 2 ) 

= log det (j N2 + p 12 H 12 H\ 2 - p^H^Ql^iHl^j 

+ log det ^I Nl + Pll H 11 H\ 1 - ^p^p^HuHlz ^/p^H u Q 12 



^PnPuH^H^ 
l P~rvQ\2^\\ 



(78) 



In 2 + Pl2Hl 2 H\ 2 yfpY2Hi2Ql2 
V / Pl2 ( 3l2-^12 Im 2 

where (a) is obtained by using Lemma 6 and Lemma 7 and (b) follows from the fact that h(Zi) = 
log det 2ire (In,), for % = 1,2. 

Fourth term: The fourth term is similar to the third term by replacing 1 and 2 in the indices. 

Fifth term: According to the fifth bound in Lemma 1, it is sufficient to replace upper bounds of 

h(Y 1 \X 2 ,Si) and h{Y 2 ) from from Lemma 7 and Lemma 5, respectively, and get the fifth bound of 
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Theorem 1 as follows 

R 1 + R 2 < h(Y 1 \S 1 ,X 2 )-h(Z 2 ) + h(Y 2 )-h(Z 1 ) 

< log det (j N2 + p 22 H 22 H\ 2 + p 12 H l2 H\ 2 + ^/ p 22 pi 2 H 22 Q l2 H\ 2 + \ / p22pi2H 12 Q\ 2 Hl 2 ^J 
+ log det 2ne (I Nl ) 

+ log det I Nl + pnHnH^ - ^/p^p^H n H\ 2 ^/pl{H u Q 12 



In 2 + P\-2.Hi 2 H 12 y/p~12~Hi 2 Qi 2 

\^Pi2Qi 2 H\ 2 Im 2 
+ logdet27re (I Nl ) - h(Z 1 ) - h(Z 2 ) 



i ,- 



y / PllPl2Hi 2 H n 
\/PllQi 2 H 11 



(b) 



log det ( I Na + p 22 E 22 E\ 2 + p 12 H l2 H\ 2 + y/p 22 pi 2 H 22 Q 12 H J 12 + v / ' p 22 p l2 H l2 Q\ 2 H\ 2 



+ log det I Nl + pnHnH^ - 



^/PiiPi2HuH{ 2 yfp^HuQ 



12 



In 2 + Pi2H\ 2 H\ 2 ^Jp~viil\ 2 Q\ 2 
Jp~v2Q\ 2 H\ 2 h 



l M 2 



-1 ,- 



^/Pllfh2H 12 H\ l 

\fp~llQ\2H\l 



(79) 



where (a) is obtained by using Lemma 7 and Lemma 5 and (b) follows from the fact that h(Zi) = 
log det 27re {In x ), for i = 1,2. 

Sixth term: The sixth term is similar to the fifth term by replacing 1 and 2 in the indices. 

Appendix B 
Proof of Achievability for Theorem 2 

In this section, we prove the achievability for Theorem 2. More precisely, we will show the following. 



Lemma 8. For a given p and H, the feedback capacity region of a two-user MIMO Gaussian IC can 



27 



achieve all rate pairs (Ri,R 2 ) G R a (H,p) such that 

R 1 < \ogdet{I Nl +p ll H ll Hl l + p 2l H 21 Hl)-N l , (80) 
R 2 < \ogdet(I N2 +p 22 H 22 Hl 2 + p 12 H 12 Hl 2 )-N 2 , (81) 
R 1 < log det (l N2 + p l2 H l2 H\^j + log det(7 JVl + Pii#ii#i f i - 

^pT 2 H n Hl 2 (I N2 + p 12 H X2 H\ 2 y x 4p-^- 2 H l2 H\^) -N 1 - N 2 , (82) 
i? 2 < log det (j Nl + p 21 H 21 Hl^j + log det(ijv a + p 22 H 22 H\ 2 - 

y/p^HnHl^In + P2iH 21 Hl 1 )- 1 ^p- 2 - 2 -p-^H 21 Hl 2 ) - Ni- N 2 , (83) 
R! + R 2 < log det (l N2 + p 22 H 22 H\ 2 + p 12 H 12 H\ 2 ^ + log det(7 JVl + PuH^H^ - 

^P^HuHI^In, + puHnHlz)- 1 Vp~^ P 1 2 H 12 h1 1 ) - - 7V 2 , (84) 
i?! + i? 2 < log det (j^ + PnHnHl + p 2 iH 2l H\^j + log det(/jv 2 + p 22 H 22 H\ 2 - 

^p 22 p 2 iH 22 Hl 1 (I Nl + p^H^H^y 1 y /p 2 ~ 2 ~p~2iH 21 Hl 2 ) -Ni- N 2 . (85) 

The rest of the section proves this result. 

Lemma 9. 77ze feedback capacity region of the two-user discrete memoryless IC includes the set of 
R 2 ) such that 

Ri < HU^X^YJ, (86) 

R 2 < nUuX^Y,), (87) 

Ri < I (Ui, Y 2 \X 2 ) + I (Xi; Y x | U u U 2 ) , (88) 

Ri < I (U 2 ; Y X \X{) + I (X 2 ; Y 2 \ U u U 2 ) , (89) 

R1 + R2 < I (Xi; Yi I £/i, U 2 ) + I (C/i, X 2 - Y 2 ) , (90) 

#1 + ^2 < / (X 2 ; Y 2 I U u U 2 ) + I (U 2 , X i; Yy) , (91) 

over all joint distributions p{ui)p{u 2 )p{x 1 \u\)p{x 2 \u 2 ). 

Proof: This result follows from Lemma 1 of [8] with U — 0. ■ 
To achieve this rate region, the authors of [8] developed an infinite-staged achievable scheme that 
employs block Markov encoding, backward decoding, and Han-Kobayashi message splitting. 

The rest of the section inner bounds this region to get the inner bound in Theorem 2. For this, we will 
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introduce some useful Lemmas. 

Lemma 10. The following holds for any Mi x Nj matrix S 

S{I Nj + 5 t 5)- 1 ^ t h 0. (92) 

Proof: It holds since it can be written as AEA^ for A = S and E = (I Nj + S^S) \ which is p.s.d. 
because E is p.s.d.. ■ 

Lemma 11. The following holds for any Mi x Nj matrix S 

det(I Nj + S^S - S^S(I Nj + S^S^S^S) < det (21 N .) . (93) 

Proof: Let us define V = S^S, we get 

det(I Nj + S f S - 5^(7^. + tfsy^S) 
= det(I Nj +S*(I Mi -SES*)S) 
= detil^+S^lM.-Sil^+S^Sy'^S) 
= det(I Nj + S j S - S^S(I Nj + S^Sf^^S) 
= det(I Nj +V- V(I Nj + Vy V) 
= det(/^. + y-y(7 iVj + V)~ 1 (V + / JV .-/ JV ,)) 

= det{I Nj +V -V{I N] -(I N] +Vy 1 )) 

= det{I N] +V{(l N] +Vy 1 )) 

= det(7 JVj + (-7^+7^+\/)((7^+\/)" 1 )) 

= det(l Nj +I Nj -(l Nj +Vy^ 

< det(27 iV ,), (94) 

where (a) follows from the the fact that V = S^S is p.s.d., and its eigenvalues are non-negative. So, the 
eigenvalues of I Nj +V are greater than or equal to 1. As a result, eigenvalues of (7jv. + V) -1 , (Ai, A^), 
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satisfy < X k < 1. So 



det^+J^-^.+V)- 1 ) = (2-Ai),...,(2-A Jv .) 



< 2^' 

= det(2/ Ar .). (95) 

This proves (93). ■ 
Now we present our achievability scheme. We use a power allocation according to 

Xi = X ip + X iu (96) 

which states that signal is divided in two parts - private, and public which are independent of each other 
with X ip ~ CN(0, K Xip ) being the private part and X iu ~ CN(0, K Xiu ) being public part such that: 

K Xip = I Mi - ^Hl(I Nj + pijHijH^y/pTjHij, (97) 

and 

K Xiu =I Mi -K Xip . (98) 

We note that the power allocation is feasible since l Ml — K x . p >z by Lemma 10 substituting ^fp~H\- 
into S. 

We will now expand the achievability in Lemma 9 using Ui = X iu for i e {1,2}. Before expanding 
each term in Lemma 9, we evaluate some entropies as follows. 

h (Yi) = log det(I Ni + puHuHl + pyjlyjl],) , (99) 

and 

h (Yi\Xi) = logdet (l Ni + P r, 1 1 ,,11,) . (100) 
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In addition, we have 



h(Yi | Ui,Uj) 
> hpilUi^Xj) 
= logdest(I N .+p ii H ii K x . p Hti 

= log det(I Ni + PuHuHl - y /p~p~H ii H\ j (I Nj + p ij H ij H^)~ 1 y/p^H ij Hi. (101) 



Moreover, we have 

(a) 

< logdet(2/ JVi ) 

= Ni, (102) 

where (a) follows from Lemma 11 by substituting ^fpJlH^ in S. This shows that h(Yi \ Uj,Xi) is upper- 
bounded by iVj. 

In our achievability, h(Yi \ Uj,Xj) appeared with a minus sign. So, without loss of generality we can 
replace it with its bound JVj for the achievability. 

The rest of the section considers the six terms in Lemma 9 and uses each of them to get the terms in 
the inner-bound of Lemma 8. 

First term: For the first term in Lemma 9, we have 

I(U 2 ,X 1 ;Y 1 ) 
= h(Y 1 )-h(Y 1 \U 2 ,X 1 ) 

( = } logdet(/ iVl +p 11 ^ii^ 1 t 1 +P2i^2i^ 2 t 1 ) -h(Y 1 \U 2 ,X 1 ) 

> \ogdet(I Nl + Pll H 11 Hl 1 +p 21 H 21 Hl 1 ) - N u (103) 

where (a) follows from (99) and (b) follows from (102). 

Second term: The second bound is similar to the first bound by replacing 1 and 2 in the indices. 
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Third term: For the third bound in Lemma 9, we have 

J(f/ 1 ;F 2 |X 2 ) + /(X 1 ;F 1 | U U U 2 ) 
= h(Y 2 \X 2 )-h(Y 2 \U 1 ,X 2 ) + h(Y 1 | U 1 ,U 2 )-h(Y 1 | U^U^X,) 

> h (Y 2 \X 2 ) - h {Y 2 \U^ X 2 ) + hiY.lU,, U 2 , X 2 ) - h {Y 1 | U 1 , U 2 , X ± ) 

= logdet (l Na + pi 2 H l2 H\^j + logdet ^I Nl + puHnH^ - p\\p\ 2 H xl E\ 2 

(I N2 + pnHuHl^yfp^HnHl^j -h(Y 2 \U 1 ,X 2 )-h(Y 1 \U 1 ,U 2 ,X 1 ) 

> logdet (^I N2 + pi 2 H 12 H\^j + logdet ^I Nl + pnHuH^ - y/pupi 2 H n Hl 2 

(I N2 + puHuHl^^pT^HuHl^j - N ± - N 2 , (104) 

where (a) is obtained from (100) and (101) and (b) follows from (102). 

Fourth term: The fourth term is similar to the third term by replacing 1 and 2 in the indices. 

Fifth term: For the fifth bound in Lemma 9, we have 

HX^Y, | U 1} U 2 ) + I(U 1} X 2 ;Y 2 ) 
= h(Y 1 \U 1 ,U 2 )-h(Y 1 \U 1 ,U 2 ,X 1 ) + h(Y 2 )-h(Y 2 \U 1 ,X 2 ) (105) 

> h(Y 1 \u 1 ,u 2 ,x 2 )-h(Y 1 \u 1 ,u 2 ,x 1 ) + h(y 2 )-h(y 2 \u 1 ,X2) (106) 

= log det(/iv 2 + P22H 22 H\ 2 + p 12 H 12 H\ 2 ) + log det ^I Nl + puHuH^ - y r p~^pl 2 ~H 11 H\ 2 

(I Na + PnHuHi^^pT^HuHi^J -h(Y 2 \U u X 2 ) -h(Y 1 \ U 1 ,U 2 ,X 1 ) (107) 

> log det(ijv 2 + p 22 H 22 H\ 2 + pi 2 H 12 H{ 2 ) + log det ^I Nl + pnHuH^ - ^ pupi 2 HuH\ 2 
(I N2 + puHnHl^^pT^Hl^j -N,- N 2 , (108) 

where (a) is obtained from (99) and (101), and (b) follows from (102). 
Sixth term: The sixth term is similar to the fifth term by replacing 1 and 2 in the indices. 
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Appendix C 
Proof of Outer Bound for Theorem 2 

In this section, we prove cross covariance matrix Q = is approximately optimal for the capacity 
region of MIMO IC with feedback. As mentioned in Section III, it is enough to prove that 



n (Q) c^o)©^,^] x [0,7V 2 ]), 

for any cross-covariance matrix Q. 

Now, we give three important relations that would be used in the main proof. 
Define E = (I N2 + ^fp~jH^(I — Q^Q\ 3 )^fp~^H\-) 1 ). The first relation is as follows 



(109) 



I Ni + puHuHl - 

\l PiiPijHiiHij \J PiiHiiQij 



— INi + PiiHii{lMi 
\fPijH- Qij 

INi + PiiHii{lMi — 



-, -1 



WjQljH 



f ^ Im< 



x / l>iil>,jlli,H,, 



Inj + PijHijHij yfp~]HijQ 

VPijQljHij 



-l _ 



/JijH, 



Q 



)H 



(a) 



IPijH\j Qi 



E —E^Jp~HijQij 
-VfhjQljHljE I M] + s/Pij Q\j HjjE^/p^Hij 



rp~ijHij 



)H\i 



(b) 



(c) 



lx + PiiHu(I - QijQlj -(I- QijQljj^fHjH^E^pTjHijil - Ch 3 (£j))ll\ 
I Ni + Pii Hu(L(I - QijQlj, v^4))4 



(d) , , 

< I^+puHiiiLil^Hj^Hl 



(110) 



where L(K, S) is as in (54), (a) follows since the inverse can be verified easily, (b) follows from finding 
the product of matrices, (c) follows from the definition of L(K, S) in (54), and (d) follows from Lemma 
4. 

The second relation is as follows 
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logdet (l N . + PijHijH^ - PijHijQijQlHl^ 
<\ogdet(l Nj +p ij H ij Hl) . 



(Ill) 



The third relation is as follows 



logdet(J JVi + puHuHl + p jill], + y/PaPjiHuQijHji + y/puPjiHjiQljHJA 

0) + + + + 

< logdet^ + puHuHl + pjiHjiH^ + paHuQuH^ + p;,;!! jiQ/uU),) 

< logdet^ + IpuHuHl + 2p ji H ji H] i ) 

< log det(/ iVj + Pn^n^ + PjiHjiHt) + iVj, 



(112) 



where (a) follows from (A - B) (At - B*) = AA^ + BB* - AB^ - BA^ y by substituting y/p^HuXt 
and yfpjiHjiXj in A and -B, respectively, (6) follows from the fact that I y Qa. 



Thus, we proved that among these three expansions, the first two expansions we started with are 
maximized by Q i3 - = while the third one is is outer-bounded by the corresponding expression with 
Qij = plus JVi. 

Now, we consider each of the six expressions in the definition of the region 1Z (Q) and outer-bound each 
expression to find the gap with 1Z o (0) being constant thus proving that 1Z (Q) C 7£ o (0)©([0, Ni] x [0, N 2 ]) 
which proves the result. 

Let the right-hand sides of the six expressions in the definition of 1Zq{Q) in (6)-(ll) be labeled as 
h(Q), h(Q), h(Q), h(Q), h(Q), and HQ) respectively. Then, the constant gap outer-bound is shown 
in the following Lemma. 

Lemma 12. We have 



HQ) < IM + N, 



(113) 



HQ) < W) + n, 



2, 



(114) 



HQ) < W), 



(115) 



HQ) < Ho), 



(116) 



HQ) < h(o) + N, 



2, 



(117) 



HQ) < HV + Ni. 



(118) 
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Proof: We start with (113). 

h(Q) = log det(J 7Vl + pnHnHl + p 21 H 21 H ] 2l + Jp^RH^QH^ + ^/p^pTiH 21 Q* H^) 

W t t s 

< logdet^ +puH ll H\ l + p 21 H 21 W 21 ) + Nj, 
= /i(0) + JVi, 

where (a) follows from (112). 

Proof of (114) is similar to (113) by replacing 1 and 2 in the indices. 
For the proof of (115) we have, 

h(Q) — l°gdet(/7v 2 + Pl2 H 12 H{ 2 - p l2 H 12 QQ^H{^j + logdet I Nl + puH^H^ - 



(119) 



In 2 + Pl2#12#12 \fP~VlH-VlQ 

\fpviQ ] H\ 2 I M2 



\J PwPyiHyiH 



ii 



VpTiQ^hI 

< logdet (l N , + Pi2#i2#i 2 ) + logdet ^ + puHuH^ - 

PuP\2H 11 H\ 2 {I N2 + P\ 2 H 12 H\ 2 )~ 1 HuHli 
= HO), 



(120) 



where (a) follows since the first expression is outer-bounded as in (111) and the outer-bound for the 
second expression can be shown on similar lines as (110). 

Proof of (116) is similar to (115) by replacing 1 and 2 in the indices. 
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For the proof of (117) we have 



i(Q) = log det (/^ + P22H22HI2 + Pi2H 12 H\ 2 + v /p22Pl2-f^22<5 t -f^l2 + y/ P22P12H12Q 

+ log det ^I Nl + puHuH^ - v /pTTpT^H 11 # 1 t 2 ^/p^H u Q ] 



In 2 + P12H12HI2 yfpviHviQ 

VpT2Q^h\ 2 h 



-. -1 



l M 2 



\J p\\p\2ii\2ii[\ 



fpiiQ^H 



11 



(a) 



< log det(ijv 2 + p22H 2 2H\ 2 + pi2H 12 Hl 2 ) + log det + pnHuH^ - 

PuPi 2 H 11 H\ 2 {I N2 + p 12 H 12 Hi 2 )- 1 H 1 2Hi^j + N 2 
= h(0) + N 2 , 

where (a) follows from (112) and using similar steps as in (110). 
Proof of (118) is similar to (117) by replacing 1 and 2 in the indices. 
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Appendix D 
Proof of Reciprocity in TZ o (0) 

In this section, we prove that replacing H and p by H and p R , respectively, and interchanging M and 
N for antennas at the nodes gives the same expressions in 1Z o (0). 
We shall prove this in two steps. In the first step we shall prove 



n (H,p) = n (H ,p R ), 



(122) 



where H = {H^, W 2l , H[ 2 , Hl 2 } and in the second step we shall prove that 



n (H',p R ) = n (H R ,p R ). 



(123) 



Clearly, the above two equalities prove the lemma. 

Let the right-hand sides of the six expressions in the definition of TZo(0) in (13)-(18) be labeled as h, 
I 2 , I 3 , 1 4, I 5 , and I 6 respectively. 
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First Step: In this step, we prove that: 



r 
1 1 


- V 




/n 


- P 

— M> 




h 


= 4 


(126) 


h 


= 4 


(127) 


h 


= 4 


(128) 


h 


= 4 


(129) 



where I' k is obtained from I k by interchanging M and JV, replacing ify with ifjj, and replacing pjj with 

Since ii and 7 3 are both bounds for R-y, I 2 and J 4 are both bounds for R 2 , and 7 5 and 7 6 are both 
bounds for Ry + R 2 , (124)-(129) will prove that TZ (H,p) = TZ {h' ,p R ). 

We start with proving (124). For simplicity we define K = (I Nl + p 2 yH 2 yH\^)~ x , K' = (I Ml + 
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P 2iH\ l H 2l )-\ and L 4 p 11 H u H\ 1 . We get 

Ji = logdet^ + pnHnHl + p 21 H 21 H j 2l ) (130) 
= log det(I Nl + p 21 H 21 H ] 2l ) + log det(I Nl + Kp lx H xl H\^) 
= logde^TT 1 ) + logdet(7 iVl + KL) 
( = } log det(K" 1 ) + log det(7 JVl + LK) 
= hgdetiK- 1 ) + logdet(/ JVl + 7.7T7) 

= logdet(^- 1 ) + logdet(7 JVl + LK(I + p 21 H 21 H ] 2l - p^H^H^)) 

= logdet(K- 1 ) + logdet(7 JVl + LK(I + p 2l H 2l H\ 1 - p 21 77 21 (7)77 2 t 1 )) 

= logdet(7<- 1 ) + logdet(7 JVl + LK{I + p 21 H 21 H j 21 - p 21 H 21 (K'^ K^H^)) 

= logde^TsT 1 ) + 

logdet(7 iVl + LK(I + p 2 i77 21 77 2 t 1 - p 2 i77 21 ((7 + p^H^H^K')!!^) 
= logdet(7<- 1 ) + 

logdet(7 iVl + LK((I + p 21 H 21 Hl) - p 21 ((7 + p 21 77 21 77 2 t 1 )77 21 K')77 2 t 1 )) 
= logdet(7<- 1 ) + logdet(7 7Vl + LK{K~ X - p 21 K~' H 21 K J )H* 21 )) 
= logdet(7 + p 21 HlH 21 ) + logdet(7 JVl + L(I - p^H^K^Hl)) 
^ logdet(7 + p 2l H\ x H 2l ) + logdet(7 + L — Lp 21 H 21 K' H^) 
= logdet(7 + p 21 77 2 t 1 77 21 ) + 

logdet(7 + pnHnHl - p 11 p 21 77 11 77 1 t 1 77 21 (7 + p^H^H^y 1 H 2 \) 
( = } logdet(7 + p 21 77 2 t 1 77 21 ) + 

logdet(7 + pnHlHn - p 11 p 2 i77 1 t 1 77 21 (7 + p^H^H^)' 1 H^Hu) (131) 
= 4 (132) 

where (a), (6) and (c) follow from Sylvester's determinant theorem [20]. (125) can be proved similarly 
due to symmetry. In addition, (126) and (127) can be obtained in the reverse direction similarly. 
We move toward the proof of (128). We should prove 

h = I' 6 , (133) 
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where 

h = log det(JAr 2 + p22H 22 Hl 2 + p 12 H 12 H\ 2 ) + 

logdet(J JVl + pnHnHl - PnPl2 H n Hl 2 (I N2 + p 12 H 12 H\ 2 )- 1 H^H^), (134) 

and 

I' 6 = \ogdet(I Ml +pnHl 1 H u + p 12 Hl 2 H 12 ) + 

log det(/ M 2 + p22H ] 22 H 22 - p 22 pi 2 Hl 2 H 12 (I Ml + p\ 2 H[ 2 H l2 y x H[ 2 H 22 ) . (135) 

If we define 



a 


_A 


log det(I N , 2 


+ P22^22H 22 


+ Pl2 H 12 Hl 2 ), 


(136) 


b 


_A 


log det(I Nl 


+ pnHnHl 


— PwP\2H\\H\ 2 (1 ' N2 + p 12 H 12 H 12 ) _1 i? 12 -f/" 11 ), 


(137) 


c 


_A 


log det(/ Ml 


+ puH^Hu 


+ p 12 H\ 2 H l2 ), 


(138) 


d 


,A 


logdet(/ M 2 


+ P22H\ 2 H 22 


— P22Pi2H\ 2 H 12 (I Ml + p 12 H\ 2 H X2 Y l H\ 2 H 22 ) 1 


(139) 



then, it is sufficient to prove a + 6 = c + d or a — d = c — b. 
Since (130) is equal to (131), we have 

logdet(7 JVl + puHuH^ + P 2\H 2 \H\^) - 

logdet(J + pnHlHn - pnp^HlMl + p^H^H^)' 1 H^H^) = 

logdet(/ + p 21 /4# 21 ). (140) 

Using similar method, we can see that 

a - d = logdet(/ Ml + p 12 H\ 2 H 12 ), (141) 

and 

c - b = log det(I N2 + p 12 H 12 Hl 2 ), (142) 

which according to Sylvester's determinant theorem [20] are equal. This proves the I 5 — I' 6 . 
(129) can be proved similar to the proof of (128) due to symmetry. 
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Second Step: It can be proved with a similar discussion as in Appendix E of [4]. A brief sketch of 
the proof is given below for completeness. 

Suppose S is a p.s.d. matrix and S* represents its complex conjugate, i.e., the matrix obtained by 
replacing all its entries by the corresponding complex conjugates. Then, it is easy to see that 

logdet(/ + S) = logdet(/ + S*). (143) 

However, note that all the terms in the different bounds of TZ o (0) are of the form of logdet(J + S). This 
in turn proves that if we replace all the channel matrices of a two-user MIMO IC with feedback by their 
complex conjugates the set of upper bounds remain the same. From this fact, it easily follows that 

n (H,p R ) = n (H R ,p R ). (144) 

Appendix E 
Proof of Theorem 5 

In this section, we will find the limit of 7£ o (0)/logSNR as SNR — > oo to get the result as in the 
statement of the Theorem 5. This follows from Theorem 3 since the capacity region is inner and outer- 
bounded by 1Z (0) with constant gaps which would vanish for the degrees of freedom. 

Before going over each of the terms in TZ o (0) and finding its high SNR limit, we first give some 
Lemmas that will be used for the proof of the Theorem. 

Lemma 13 ( [4]). Let e C NjXMi be a full rank channel matrix. Then, the following holds 

log det (l N . + PijHijH}^ = aij mm (M h Nj) log SNR + o(log SNR ). (145) 



Lemma 14. Let E e C NxM be a diagonal matrix with elements a±, ...,a m where m = min(M,N) and 
A G C mxm be a diagonal matrix with elements |o"i| 2 , \cr m \ 2 , then 



C^ + A)- 1 

I(N-M)+ 

I m -(I m + A)- 1 

O(M-N) 



(146) 
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Proof: We will split the proof in two cases, depending on whether M > N or M < N. 
Case 1 - M > N: In this case, we have 



(4* + A) -1 





a\ 

'•• 





kil 2 
i+kiP 











i+kil 2 













l (N-M) + 



1 



i+kil 2 













1 + Wm\ 2 



a i 
'•• 

a m 






Wm\ 2 
l + km| 2 









0(M-A0+ 





1 - 






1 

i+k m | 2 








0(M-JV)+ 



/ m -(/ m + A) 1 





0(M-JV)+ 



(147) 
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Case 2 - M < N: In this case, we have 



(J m + A)- 1 
/. 



(N-M) + 
1 



a{ 
'•• 


kil 2 
i+kiP 








a* 






1 

l + |<Tl| 2 










\<Tm\ 2 
1+kml 2 





i+kil 2 




















1 

i+kml 2 








T(N-My 

















o m 











1 - 




0(M-AT)H 












i+kml 2 

0(M-AT)H 



I m -{I m + X)- 1 

0(M-Ny 



(148) 



Lemma 15 ( [4]). Let Ha e C NiXMi and Hji e £ NtXMj be two full rank channel matrices such that 
[HaHji] is also full rank. Then, the following holds 

fogdestiln+piiHiiHl + pjiHjiHt) = f(N t , (a ii: M,) , (a ji: M,-))logSNR +o(logSNR) (149) 

where f is defined in 28. 

Lemma 16. Let Ha e C NiXMi and Hij e <£ N i* M j fr e p^ Q channel matrices with each entry independently 
chosen from CN(0, 1). Then, the following holds with probability 1 {over the randomness of channel 
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matrices). 



log det{I Ni + puHaHl - y/PiiPijHiiH^Ijf. + />;, I I , , / /, ; ) v V„7^7/, ; 7/„) 

= [c^min ((Mj - A^) + , ty) + («« - a tj ) + (min (M,, A^) - min ((M, - A^) + , N t ) )] logSNR 

+o(logSNR). (150) 

Proof: Let the singular value decomposition (SVD) of the channel matrix be given by = Vij 
^ijUjj, where G U NjXNj and CT^- e JJ MiXMi are unitary matrices and e U NjXMl is a rectangular 
matrix containing the singular values along its diagonal. Using the SVD of the matrix 77^ we get 



Ay. + PiiHuH^ - y /p~p~H ii Hj j (I Nj + p. Jljjll- A 1 v /p ii p ij H ij H ii 

I Ni + pa II;; ^I Mi ~ PijH-j (jNj + PijHijH}^ llj^j II ;; 



(151) 



(a) 



1^ + PuHh(Im i — PijHijVij 



7jV; + Pit Hu ( I Mi — PijUijUlj 



(J^.+SNR^A^)- 1 

^(Nj—Mi) 

(7 roi .+SNR a «A ij )- 1 

/, 



= 7^ + p ll H ll U l] {I Mi - SNR Q -4 



(6) 



7jv 4 + PiiHaUij^Mi 



(7^+SNR^A^)- 1 
/„,.• (/„, • SNR" A; ; ) 1 











= I Ni + SHR aii H u U, 



(7 roi .+SNR a «A iJ -)" 1 







(152) 



where (a) results from SVD of the matrix 77^ and (b) follows from Lemma 14. 

Let us decompose e JJ MiXMx into two parts, C/yi and C/jj 2 such that C/y = [C/yi C/^], where 

Uij! G £/^xmin{M^ 3 } and ^ e jjMi x (M i -iV J )+ _ w£ get 
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logdet^ + SUR aii HuiUij 



[ (M i -A r j ) H 



logdet(/ A r i + PiiHuHly/piiPijHuHlj (I Nj + p, s 1 1 , ; / / • ; ) 1 y/p ii p ij H ij Hl i ) 

(J^.+SNR^A^)- 1 
J ( . 

= logdet(7 JVj + /^(SNR^C/^/^. + SN^A^) -1 ^ + SNR^E/^E^iZl) 

= logdet(/ JVi + SNR^tf^C/^tft + SNR^-^^^^SNR-^/^. + K^U^Hl (153) 

where = min(Mj, A^), A i;? - is a diagonal matrix containing the non-zero eigenvalues of H i3 B.\ 3 . 

We note that A^- is invertible and when SNR is large, we can bound SHR~ aij I m .. + A i3 from above and 
below as, A i3 - ^ SNR _Q!lJ : I mij + A„ ■< I + Ay. We will only pursue the direction where SNR~ aiJ 7 mij + 
Aj-, >: h-ij and can see that both the directions produce the same result and thus replacing the inner 
and outer bound by equality. In what follows, even though SNR~ aiJ I m + A i3 >z A i3 , we will substitute 
SNR~" iJ ' I mij +Aij = Aij since by the inner and outer-bounding approach, it can be seen the the limit will 
be exactly the same thus not causing any difference in the result. Thus, we have 



logdet(/ 7Vi + p,, H n ///, sjpapij II a ( /.v ; + Pi 3 Hi 3 H] 3 ) 1 y /p ii p ij H ij Hl i ) 
= logdet(7 A r i + SHR^H^Uj^Hi + SUR^~ a ^ H^iK^U}^ + «(logSNR) 
( = } logdet(7 A r i + SHR aii H u U ij2 Uj j2 Hl + Sm^~ a ^ + H^U^A^U^Hl 
( = } /(^,(« ll ,(M l -iV J ) + ),(K-a y ) + ,mm(M J ,iV J )))logSNR +o(logSNR) 
= [c^min ((M, - N 3 ) + , + [a u - a l3 ) + min((A^ - (M, - A^) + ) + , N 3 , M 4 )logSNR 

+o(logSNR ) 

( = } [c^min ((M< - N 3 ) + , N,) + (a lt - a l3 ) + (min (M i: N t ) - min ((M< - iV j ) + , A^) )] log SNR 
+o(logSNR ), (154) 

where (a) follows from the fact that if (an — a i3 ) is less than zero we have 

Sm^-^HuU^A^U^Hl = o(logSNR), (155) 

(b) follows from Lemma 15 and that HaU^i, HaUijiA' 1 ^ 2 and H u \U i32 U i3 iA~^ 2 \ are all full rank with 
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probability 1; (c) follows from some simple manipulations. ■ 
The rest of the section considers the 6 terms in TZ o (0) in (13)-(18), and finds the GDoF region for the 
MIMO IC with feedback. 

First term: According to the first bound in 1Z o (0), we have 

logdet(/ JVl + pnHuHl-L + p 2X E 2l li\^) 
( = } /((iV 1 ,(a 11 ,M 1 ),(a 21 ,M 2 )))logSNR + o(loglogSNR ), (156) 

where (a) is obtained from (15). Now, dividing both sides by logSNR, we get the first GDoF expression. 
Second term: The second bound is similar to the first bound by replacing 1 and 2 in the indices. 
Third term: According to the third bound in 1Z o (0), we have 

logdet (l Na + p 12 H 12 H\ 2 ^ + logdet(/ iVl + puHuH^ - 
\Zp~yvPviH 1 iH\ 2 (I N2 + pnH^H^y 1 ^/p^p^HuH^) 
= a 12 mm (Mi, N 2 ) + cmmin ((Mi - iV 2 ) + , JVi) + 

(an -tti 2 ) + min(Mi,A^ 1 ) - min ((Mi - iV 2 ) + , N ± ) +o(logSNR), (157) 

where (a) is obtained from Lemma 13 and Lemma 16. Now, dividing both sides by logSNR, the third 
GDoF bound results. 

Fourth term: The fourth term is similar to the third term by replacing 1 and 2 in the indices. 
Fifth term: According to the fifth bound in R o (0), we have 

log det (/jv 2 + p22H 22 H\ 2 + p 12 H 12 H\ 2 ^j 

+ logdet(i A r 1 + pnHuH^ - ^pnpi2H u Hl 2 (I N2 + p 12 H 12 H\ 2 y l ^ ' p\\p\ 2 E\ 2 E[ x ) 
= /((iV 2 ,(a 22 ,M 2 ),(ai 2 ,Mi))) + ai 1 min((Mi-iV 2 ) + ,iV 1 ) + 

(an -a i2 ) + (min(Mi,iVi) - min ((Mi - N 2 ) + , iV\) ) + o(logSNR ), (158) 

where (a) is obtained from Lemma 15 and Lemma 16. Now, dividing both sides by logSNR, the fifth 
GDoF bound results. 

Sixth term: The sixth term is similar to the fifth term by replacing 1 and 2 in the indices. 
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